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ABSTRACT 

Between 1994 July and 1997 February we monitored the optical spectrum 
of RX J0019. 8+2156. This supersoft X-ray source is one of only two accreting 
white dwarfs in the Galaxy that are thought to be burning hydrogen on their 
surface as a consequence of a high rate of mass transfer from a binary companion. 
Accurate orbital ephemerides are derived from radial velocity measurements 
for the white dwarf, which are obtained from strong Hell emission lines with 
a stable velocity semi-amplitude oi K — 71.2 ± 3.6 km s~^. We report the 
discovery of transient, low-velocity, bipolar jets. These jets are respresented by 
redshifted-blueshifted pairs of emission lines from H and Hell with an outflow 
velocity of v cos(i) ~ 815 km s~^, where i is the binary inclination angle. When 
present, the jet lines seen in Ha also exhibit an orbital modulation of 71 km 
s~^ which strengthens the interpretation that this is the orbital velocity of the 
white dwarf and also indicates that the jets are oriented nearly perpendicular to 
the orbital plane. On most occasions, the H emission line profiles are further 
altered by P Cygni absorption effects, and the strength of this absorption is also 
dependent on binary phase. We show that the jets and the P Cygni features 
have very different temporal characteristics and binary phase dependence; thus, 
we conclude that the outflowing material and the absorbing wind must have 
essentially different geometries. Finally, the measured mass function is combined 
with binary evolution models to suggest a limit on the inclination angle, i < 40°. 
A particular model invoked to explain a high rate of mass transfer requires 
16° < i < 25°. However, at such small inchnation it is difficult to explain the 
large amplitude of the orbital light curve (~ 0.5 mag). Alternatively, the results 
could signify substantial vertical structure in the accretions disks of supersoft 
X-ray sources. By contrast, a simple model fit to the jet outflow lines indicates 
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an orbital inclination angle of 35° < i < 60°. 



Subject headings: stars:binaries:close, X-rays:stars, starsrcataclysmic variables 
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1. Introduction 

Supersoft X-ray sources exhibit extremely steep spectra in the range of 0.1 - 2 keV. 
They can be divided into three general subclasses: low and high luminosity galactic objects, 
and AGN. In this paper we concern ourselves with luminous galactic supersoft X-ray sources 
(hereafter SXS), which have characteristic luminosities of ~10^^ - 10^^ erg s~^ and effective 
temperatures of ~2-6xl0^ K {kT ~ 17 — 50 eV). The first luminous supersoft sources were 
discovered during a survey of the Large Magellanic Cloud with the Einstein satellite ( [Long 



et at. 195I| ). The ROSAT X-ray survey ( |Voges 19931 ) has now established that the SXS 



constitute a distinct astronomical class. Greiner (1996) has cataloged and summarized the 
properties of 34 known SXS. Only seven of these objects are in the Galaxy, while seven are 
in the LMC, four are in the SMC, and sixteen have been located in M31. 

A particular model for at least some of the long-lived SXS invokes steady or cyclical 
nuclear burning of matter on the surface of a ~1 Mq white dwarf that is accreting mass from 
a main-sequence or subgiant companion star of ~ 1 — 2 Mq at a rate of 3x10^^ — 10^^ Mq 
yr^^ ( [van den Heuvel et al. 19"92| , hereafter HBNR; [Rappaport, DiStefano, fc Smith 1994 



hereafter RDS; and references therein). Such high transfer rates are a natural consequence 
of thermal-timescale mass transfer via Roche lobe overflow, which characteristically occurs 
in this type of system. The transfer is potentially unstable because the donor star is more 
massive than the accreting white dwarf; however, the rate of transfer is limited by the 
thermal time scale of the donor star. 

Models for the evolution of these systems (RDS) indicate that the orbital periods 
should lie in the range of 8 hr to 1.4 d. Moreover, these studies predict that there should 
be some 10^ such systems in the Galaxy and in M31. In spite of the greater distance, it 
is apparently easier to detect these supersoft sources (at least in the soft X-ray band) in 
nearby external galaxies at relatively high 6^^, such as the LMC (—33°) and M31 (—21°). 
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In the plane of our own Galaxy, a hydrogen column density of only ~10^^ cm~^ will absorb 
most of the very soft X-radiation. In the Galactic plane this corresponds to distances of less 
than a kpc. 

Five of the seven SXS in the Milky Way are ttnlikely to be the type of nuclear burning 
white dwarf (NBWD) discussed above. Two are recently exploded novae, two are symbiotic 
novae, and one is unclassified. Thus, there are only two good NBWD candidates in the 
Galaxy: RX J0925.7-4758 QMotch, Hasinger, fc Pietsch 1994|) and RX J0019.8+2156 



(hereafter RXJ0019; [Beuermann et al. 1995|) . Since RX J0925.7-4758 is obscured by the 



Vela Sheet molecular cloud, and its intrinsic optical spectrum is heavily reddened ( [Motch, 



Hasinger, fc Pietsch 1994 ; Grabelsky et al. 1987| ), RXJ0019 remains as the one prime 



galactic NBWD candidate for detailed optical study. 

Beuermann et al. (1995) discovered RXJ0019 as part of the ROSAT all-sky survey in 
1990 December. They subsequently initiated a battery of short ROSAT PSPC, ROSAT 
HRI, lUE, and optical observations in 1992-93 ( [Beuermann et al. 1995| ). The following 



is a brief summary of their results. The optical counterpart to RXJ0019 is a moderately 
bright object with V = 12.2. The light curves of the X-ray, optical, and UV data show 
modulations with a period of 15.85 hr, interpreted as the binary period, while smaller, 
aperiodic fluctuations are seen on time scales of hours. They adopt a distance to RXJ0019 of 
2 kpc on the basis of the calcium absorption features and the X-ray column density {Nh ~ 
4.4x10^° cm~^). The X-ray spectrum is extremely soft and luminous. For an assumed LTE 
model atmosphere ([van Teeseling, Heise, fc Paerels 1994|), the ROSAT spectrum yields 



kT = 34 eV and an estimated bolometric luminosity, L^oi — 4x10^^ erg s~^ with significant 
uncertainty caused by the limited spectral constraints. The optical spectrum shows strong 
hydrogen lines, exceptionally strong He II emission, and P Cygni absorption associated with 
the Balmer lines. Beuermann et al. also report that the radial velocity of the He II A4686A 
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line is modulated at the orbital period with a semi-amplitude of 67 ± 3 km s . 

Greiner & Wenzel (1995) examined archival photographic plates dating between 1890 
and 1990 from the Harvard plate stacks and the Sonnenberg Observatory. They detected 
the orbital period, low amplitude variations on time scales as short as 2 hr, and ~1 mag 
variations over intervals of ~20 years. 

Will & Earwig (1996) studied RXJ0019 with UBVRI photometry from 1992 to 
1995 using a multichannel photometer on the Wendelstein Observatory 80-cm telescope. 
RXJ0019 shows both primary and secondary minima; the primary minimum is a ~ 0.5™ 
eclipse which covers nearly 40% of the orbital phase (FWHM). Will & Earwig's analysis of 
mid-eclipse times provides an accurate orbital ephemeris with a period of 15.85133 ±0.00017 
hr. 

Gansicke, Beuermann, & deMartino (1996) observed RXJ0019 with the lUE telescope 
using both the short-wavelength prime (SWP; 1250-1900 A) and long wavelength prime 
(LWP; 2000-3200 A) spectrographs. The 15.85 hr periodicity, and both primary and 
secondary dips are present in the UV light curves; primary minima were deeper in the LWP 
spectra (32% modulation) than in the SWP spectra (22%). The Hell 1640A line is present 
and shows a stronger modulation than the continuum (43%). Fits to the UV spectrum 
using a simple model for an irradiated donor and accretion disk ( PR^einsch et al. 1996a| ) yield 
parameters consistent with masses of ~1.0 and ~1.5 M© for the white dwarf and donor 
star, respectively QGansicke, Beuermann fc deMartino 1996| ). 

In this work we present our spectroscopic observations of RXJ0019. The radial velocity 
work of Beuermann et al. (1995), was completed in four days. In contrast, our observations 
span ~ 2.5 yr and allow us to probe the spectral features and the stability of radial velocity 
measurements over a much longer baseline. In §2, we present our observations and the data 
reduction. In §3 we present the results from the spectral analysis and the radial velocity 
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study, including the P Cygni profiles and transient high-velocity emission lines, which we 
interpret as emission from bipolar jets. Similar high-velocity lines have been seen in recent 
infrared spectra reported by Quaintrell & Fender (1998). In §4 we discuss our attempts to 
model the complete binary system using the constraints of these and previously published 
data on RXJ0019. 

2. Observations & Data Reduction 

We monitored RXJ0019 at the Michigan-Dartmouth-M.I.T. (MDM) observatory on 
Kitt Peak between 1994 December and 1997 February. Most observations were made with 
the 1.3™ telescope, the Mark III spectrograph, a 600 //mm grism, and a CCD detector. In 
most cases, our instrument configuration yielded a pixel size of 1.4 A and a resolution of 
~3A (FWHM). On occasion we used a different CCD with wider pixels and a resolution 
(FWHM) of ~5A (see Table |l|). The 1995 November observing run was made on the 2.4™ 
telescope with similar instrumentation and ~3A resolution (FWHM). In all, we acquired 
757 spectra of RXJ0019 with a total exposure of 180 ks. Table |I] summarizes the MDM 
observations. 

We also monitored the source at the Fred Lawrence Whipple Observatory (FLWO) 
located on Mt. Hopkins, AZ between 1994 July and 1995 February. We used the 1.5™ 
telescope and the FAST spectrograph ( [Fabricant fc Cheimets 1994| ) and a pair of 1200 //mm 
gratings. The twenty-five FLWO observations, which are listed in Table ^, were of longer 
duration and of higher spectral resolution (~2A FWHM) than the MDM observations. To 
maintain a moderate range of wavelength coverage, these data were obtained in alternating 
portions, one red (5750-6700 A) and one blue (4550-5500 A). 

Given that the orbital period of RXJ0019 is 15.85 hr, it is not possible to observe one 
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complete cycle in a night. Furthermore, as the period of RXJ0019 is in a near 2:3 resonance 
with the day/night cycle, it is difficult to obtain uniform phase coverage. However, these 
shortcomings are partially circumvented by our 2.5 year observing campaign. During a 
typical MDM observing run, we pointed at RXJ0019 several times per night (roughly every 
two hours) , and we obtained five 200-s exposures per pointing. Occasionally we continued to 
take 200-s exposures of the target for several hours to search for effects such as non-orbital 
fluctuations in the velocities of the emission lines. 

We used standard IRAF routines to reduce the raw CCD images into one-dimensional 
extracted spectra and also to calibrate the spectra in wavelength and flux. In performing 
the wavelength calibrations, we used the lamp spectra taken immediately before and after 
each pointing sequence to interpolate the dispersion solution to the mid-exposure time of a 
given observation. 

3. Analysis and Results 

3.1. Radial Velocity and Period Determination 

Radial velocity measurements of RXJ0019 were extracted with IRAF tools that 
compute cross correlations between two continuum-subtracted spectra. In practice, the flrst 
of these spectra is a "template" , which is selected to have very nearly the same spectral 
features as the object spectra. The velocity of an object spectrum, relative to a restframe 
which is deflned by the template, is determined by cross-correlating the two spectra. All 
of the MDM observations of RXJ0019 (1.4A/pixel) were cross-correlated with a template 
spectrum constructed from the average of the 1996 January observations. Because there 
is P Cygni absorption associated with the Balmer lines (see below), the cross correlations 
were limited to wavelength intervals selected for Hell hues located at AA4200, 4542, 4686, 
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5412, and the series in the range of 6004—6407 A. A similar analysis was conducted with 
the FLWO data, using the mean FLWO spectrum as the cross correlation template. Finally, 
prior to radial velocity determinations, each of the two template spectra was artificially 
redshifted so that the He II line centers were maximally consistent with their corresponding 
laboratory wavelengths. The Hell line profiles in RXJ0019 are relatively narrow and 
symmetric, and we do not expect significant systematic errors to arise from this procedure. 

For the MDM spectra, the radial velocities were computed and then screened for 
discrepant values by flagging results that are statistically inconsistent (> 3a) from the 
other 4 measurements obtained during the same telescope pointing. Only 19 spectra were 
eliminated for this reason. Additionally, one set of observations (December 1995b), produced 
unusually large variances within groups of velocity measurements, and we eliminated these 
60 spectra from further analyses. Thus, we derived a total of 691 velocity measurements 
from the two observatories. We used the mid-exposure time of each observation to convert 
the velocities to the heliocentric rest frame before proceeding with the timing investigations. 

The binary period search was performed with a Stellingwerf dispersion minimization 
routine (^tellingwerf 1978| ). The variance statistic, B, is essentially the ratio of the average 



variance computed in small phase bins to the variance of the data as a whole. When the 
data are randomly distributed, G is approximately unity; however, if the data become 
ordered, as when folded at the orbital period, G is reduced significantly. 

The 691 velocity measurements were included in the search for the period that 
minimizes G. The search covered trial periods in the range of 8 - 20 hr. The results are 
shown in Figure |^, and the ~15.85 hr period is easily identifiable. The other two highly 
significant features are related to the orbital frequency (~1.514 cycles per day): the dip at 
9.5 hr corresponds to the 1-day alias period (or 2.514 cycles per day), while the dip at 19.1 
hr corresponds to a frequency that is one-half of the one-day alias. 
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We computed a least-squares parabolic fit for the values of to localize the minimum 
near 15.85 hr. The results are shown in Figure]^, and the derived spectroscopic period and 
uncertainty are included in Table |^. We folded the radial velocities about the spectroscopic 
period, as shown in Figure ^ and then we fit the results to a sine function plus a constant 
to determine the best values for the velocity semi-amplitude, the systemic velocity for 
the binary system, and the time of maximum velocity. The best values for the binary 
parameters are given in Table |^. 

Our spectroscopic period is consistent with the comparably accurate photometric 
period determined by Will & Barwig (1996; see Table ^). However, there is a discrepancy 
regarding the relative phases of the velocity and photometric modulations. According to 
Will & Barwig (1996), a primary minimum should occur at HJD 2450014.273 ± 0.007, 
which corresponds to a spectroscopic phase of 0.171 ± 0.013. However, if one assumes that 
the binary orbit is circular and that the accretion disk is circularly symmetric, then the 
photometric minimum should occur at the time of superior conjunction of the accreting 
star. This conjunction corresponds to a spectroscopic phase of 0.25. The phase offset 
between spectroscopic and photometric results is therefore 0.079 ±0.013, a discrepancy with 
a significance of 6cr. Using data from 1992 August, Beuermann et al. (1995) also reported 
that photometric minimum occurs 0.07 in phase earlier than one would expect from the 
spectroscopic ephemeris, so this phase offset appears to be a stable characteristic of the 
binary system. 

3.2. Deep Optical Spectrum 

Our extended monitoring provides not only radial velocity measurements, but also 
about 160 ks of selected exposure at a resolution of ~3A. While peculiarities in the spectrum 
of RXJ0019, e.g., the extremely large He II(A4686)/H/3 ratio, have been known since its 



- 11 - 



discovery, the high signal-to-noise of this new data set allows us to discern previously 
unreported spectral features. 

The different instrumental configurations have different spectral bandpasses, but a 
substantial subset of the MDM observations has a common region of overlap in the range 
of 4280-7000A (see Table p. We constructed a "deep rest frame" spectrum using the 398 
observations obtained during 1995 June, 1995 December, 1996 January, and 1996 December. 
Each of the flux-calibrated spectra was corrected for the Doppler motion predicted by the 
spectroscopic ephemerides (Table and the average of these results is the rest frame 
spectrum shown in the top panel of Figure ^. The P Cygni absorption features at H/5 and 
H7 are clearly evident, as is the strength of the Hell lines (e.g.. Hell A4686) relative to 
lines of He I or H. 

To enhance the visibility of subtle features in the rest frame spectrum, we fitted the 
continuum to a smooth function and then subtracted the continuum from the rest frame flux 
densities. The results are shown in the bottom panel of Figure ^ for the wavelength interval 
of 5600-6450A. This spectrum shows the n-^b transition series in He II, while the stronger 
He II lines in Figure ^ arise from n— i>3 and transition series. To our knowledge, this is 
the most detailed set of Hell features observed in any celestial object. In Table ^ we list all 
of the emission lines detected in the rest frame spectrum, along with the line identifications 
(when known) and the measured equivalent widths. The high ionization levels represented 
by these spectral features are very likely caused by photoionization related to the luminous, 
supersoft X-ray emission in this system. 
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3.3. Transient Jets 

In addition to the strong emission lines and P Cygni Balmer profiles, some of the 
spectra of RXJ0019 exhibit red-shifted and blue-shifted pairs of high-velocity lines. One 
such example, shown in Figure ||, is the red and blue sequence of 900 s exposures obtained 
at FLWO on 1994 December 12. The high velocity lines are clearly separated from the 
central line, and we interpret these results as the discovery of transient episodes of bipolar 
outflow in RXJ0019. High- velocity line pairs are evident at Ha, H/?, and Hell A4686 with 
approximate equivalent widths of 3.0, 0.5, I.OA, respectively. These lines are transient 
on time scales of months. We also observed them during the MDM observing runs of 
1994 December, 1996 May, and 1996 December, while they were not detected during the 
remaining observing runs. RXJ0019 is now the second supersoft X-ray source to show 
high-velocity emission lines. The LMC source RX J0513. 9-6951 shows weak pairs of emission 
lines near Hell A4686 and H/5 corresponding to an outflow velocity of approximately of 
3900 km s"^ ( |Crampton et al. 1996| ; [Southwell et al. 1996|) . 

To track the radial velocities of all the spectral components related to Ha, we first 
subtracted the continuum from the spectra within ±4000 km s~^ of the line position. 
We then applied a least-squares fit to the line profiles, using Gaussian functions to 
represent, whenever present, the central line, P Cygni trough, red outflow, and blue outflow 
components. In Figure |^ we plot the centroids of the Gaussian fits to the central and 
high- velocity lines as functions of the orbital phase for Ha. The radial velocities of these 
three components move in unison. The results of chi-squared fits to these curves are given 
in Table ^. All three Ha lines are in phase with the He II velocity curve to within statistical 
uncertainties. Furthermore, we find that the K amplitudes of the central, red, and blue 
lines of Ha agree with the i^-velocity of the Hell A4686 line to within 2.0o", O.lcr, and 1.6cr, 
respectively (see Table |). 
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The presence of P Cygni absorption (and any phase dependence in its strength) will 
introduce systematic errors into the systemic velocity (and fC-velocity) values derived for 
the central component and the blue outflow line. Both of these effects are seen in Table ^ 
On the other hand, the red outflow line, unaffected by P Cygni absorption, exhibits a 
velocity curve in excellent agreement with the Hell results. After correction for the 
systemic velocity of -57.5 km s~^ determined from the Hell lines, the red outflow line has 
a mean velocity of 815.5 ± 10.0 km s~^, while the blue outflow line has a mean velocity 
of —851.9 ± 14.9 km s~^, which is predictably shifted downward by P Cygni absorption 
associated with the central component. We conclude that the jet lines are consistent with a 
symmetric model in which the outflow velocity is 815 km s^^. 

We note that the central and outflow lines at Ha are all significantly broader than the 
spectral resolution of the FLWO instrument (~90 km s~^ at Ha). We measure an intrinsic 
width (FWHM) of about 400 km s~^ for the jet lines. This information is used below to 
constrain geometric properties of the bipolar outflow in RXJ0019. 

Finally, in Figure we plot the equivalent width of P Cygni absorption versus orbital 
phase. There is a clear variation in the strength of the absorption feature as the binary 
system rotates. In particular, there is an absence of P Cygni absorption from spectroscopic 
phase 0.8 to about 1.0, which corresponds to the phase quadrant that begins just after 
the white dwarf passes closest to the observer and ends when the two stars are aligned 
perpendicular to the line of sight. The orbital phase dependence accounts for much of the 
hour-to-hour variation in the P Cygni profile noted above. However, the scatter in Figure |^ 
indicates that there is still a significant random component to the absorption column. We 
return to a discussion and interpretation of the P Cygni profiles in §4. 
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4. System Modeling 

In this section, we focus on four specific results of the optical spectral analysis to 
deduce additional binary system parameters and to construct a simple kinematic model for 
the outflow of material from the vicinity of the accretion disk. The spectral results that we 
utilize are summarized below: 

1. The Doppler velocity curves for Hell (A4686), Ha, and the red and blue jet lines of 
Ha are in phase and have a common K amplitude of 71 km s^^. One expects the high 
excitation He II lines, the central Ha line, and the jet lines to all originate in different 
physical regions. The most reasonable interpretation for their common i^-velocity is 
that it reflects the orbital motion of the compact object. 

The literature on cataclysmic variables contains many warnings that the velocity 
curves derived from H emission lines do not yield reliable dynamical results (e.g., 
[Shafter, Szkody, fc Thorstensen 1986| ). In such cases the line proflles indicate a 



complexity of emission-line sources, and the derived K velocities depend on the 
manner in which the cross correlations are computed. Nevertheless, in RXJ0019 there 
are strong arguments in favor of the interpretation that the measured i^-velocity does 
represent the orbital motion of the white dwarf: First, our exclusive use of Hell lines 
minimizes contamination by low-excitation emission components far from the white 
dwarf. Secondly, we see no evidence of systematic changes in the K velocity over a 
2.5 year time scale; moreover, our results are also consistent with the earlier work of 
Beuermann et al. (1995). Finally, as noted above, we would not expect the central 
and outflow velocities to exhibit the same velocity modulations unless there were a 
common origin for their motion. Our interpretation should, of course, be reassessed 
as future optical measurements are made available. 
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2. The red and blue satellite lines, when present, have mean projected velocities of 815 
km s~^ in opposite directions. 

3. The intrinsic width of these outflow lines is ~400 km s~^. 

4. P Cygni profiles are associated only with the Balmer lines, and are most pronounced 
on the H/? and higher transitions. The strength of the P Cygni profiles varies 
systematically with orbital phase. 



4.1. Monte Carlo Constraints on Binary System Parameters 

We utilize a model-dependent Monte Carlo method for estimating a number of the 
binary system parameters for RXJ0019, starting with only the known orbital period (P) 
and the K amplitude (71.2 ±3.6 km s~^) of the accreting compact object, which is presumed 
to be a white dwarf. The mass function can then be expressed as: 



Mdon or sin^ i ^ K^P 
Mmd\'^ 2ttG 



f{m) = — 2 = ^177 = 0-0247 ±O.OO4OM0. (1^ 

(l + ^rf^' 

\ donor / 



It is apparent from the small value of the mass function that in order for the donor star 
to have a plausible mass, either the inclination angle must be small, or the mass of the 
accretor must be large compared to that of the donor. 

If we utilize a few model-dependent ideas about how an accreting binary such as 
RXJ0019 may evolve, we can combine the measured mass function with a Monte Carlo 
sampling of evolving binaries to infer a substantial number of system parameters. Such an 
approach was developed previously to infer the system parameters for the binary X-ray 
pulsar GROJ1744-28 QRappaport fc Joss 1997|) . We briefly outhne that method here. 



We start with a short period binary system consisting of a donor star and a white 
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dwarf with initial masses in the following ranges: 



0.8 < Mdonor/MQ <3.0 



(2) 



0.6 < M^d/MQ < 1.2 . 



(3) 



We begin the binary evolution calculation only after the formation of the white dwarf. 
This occurs at the end point of a common envelope phase in which the giant progenitor of 
the white dwarf has its envelope stripped off by the companion star (the "donor" in the 
supersoft X-ray source phase). The lower limit of ~0.8 M0 for the mass of the donor star 
is based on the requirement that the binary system eventually reaches a period of 15.85 
hr. Donor stars with lower initial mass will come into Roche-lobe contact while on the 
main sequence with binary periods much shorter than 15 hr and will evolve toward shorter 
orbital periods (except in the case of extreme X-ray heating effects; see discussion below). 
In initially wider orbits, such low-mass donors would have insufficient time within the age 
of the Galaxy to evolve and fill their Roche lobe. The lower limit for the mass of the white 
dwarf (0.6 M0) is required to attain the high surface temperatures observed in supersoft 
X-ray sources ( [Iben 198^ ; HBNR; RDS). Finally, the upper limit of 1.2 Mq for the mass of 
the white dwarf is taken from the population synthesis study of RDS. 

Our approach is to choose, via a Monte Carlo method, the post-common envelope 
properties of the binary system, evolve the binary to the point where its orbital period 
reaches 15.85 hr, and then find the orbital inclination angle that will yield the observed K 
velocity of the white dwarf. If the orbital period never reaches 15.85 hr, then the system 
is discarded. If the system does reach a state in which it could be a plausible match for 
RXJ0019 (i.e., with a 15.85 hr period and an orbital velocity for the white dwarf at least 
as large as the observed i^- velocity), then the system parameters at that point are saved. 
The cumulative statistics for all the binaries that produce RXJ0019-like systems form 
probability distributions for the orbital inclination angle, donor mass, core mass of the 
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donor, and the luminosity of the donor. We may then impose further constraints on the 
RXJ0019-hke systems, such as the mass transfer rate. This yields probability distributions 
for the system parameters that are consistent with the mass transfer rate being sufficiently 
high to support steady (or quasi-steady) nuclear burning on the white dwarf, ~ 3xl0~® 
< M < 10-6 Mq yr-i ( |lben 1982| ; HBNR; RDS). 

In principle, the best way to choose the initial {i.e., immediately after the common 
envelope phase) binary system parameters is from the output of a population synthesis 
code (see, e.g., RDS). However, for purposes of this exercise it is sufficient to choose the 
masses of the white dwarf and the donor star with uniform probability over the ranges 
that are permitted. We must also specify the initial core mass, Mc, of the donor star. For 
simplicity, we pick the core mass from a uniform distribution over the range < < 0.15 
Mq, above which point the orbital period of the system would always be longer than 15.85 
hr. The evolution is started with the stars separated by a distance which corresponds to 
the secondary star just filling its Roche lobe. Further details about the binary evolution 
code are given by DiStefano et al. (1997). 

This procedure of choosing initial binary parameters and then evolving the system 
to see if it would become a near match to RXJ0019 was carried out 2x10^ times. We 
recorded the binary parameters for 437,000 successful systems. In each case the results were 
weighted by a factor of sin^i/ cosi, which takes into account the probability of observing 
randomly oriented systems with the correct inclination angle to produce the observed mass 
function (see [Rappaport fc Joss 1997| ). The results are shown in Figure ^ where probability 
histograms are given for the inclination angle and mass of the donor star. The current 
value of the binary period of RXJ0019 can be approached from either longer periods (with 
P < 0), or after the binary has passed through its minimum orbital period (with P > 0), 
and we keep track of these cases separately. Each histogram in Figure ^ also displays 
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shaded regions that indicate those binaries with moderately high rates of mass transfer 
(3x10^^ M0 yr~^ < M < 10^^ Mq yr^^; single hatching) and those with very high rates of 
mass transfer (M > lO^^Mo yr~^; double hatching). The unhatched histograms include all 
mass transfer rates. 

From this figure we learn that the values of binary period and the emission-line 
K- velocity, combined with acceptable values for the initial stellar masses, constrain the 
inclination angle to the range of 16° — 40°, regardless of the mass transfer rate. Furthermore, 
if we assume that there is at least a modest rate of mass transfer (M > 3x10^*^ yr~-^), 
then the inclination angle is nnlikely to be greater than 30°. Both of these inclination limits 
can be avoided only if (i) the K velocity that we have used is not the correct one to be 
associated with the orbital motion of the accreting star, or (ii) one allows for more extreme 
constituent masses, e.g., a 1.5 Mq donor star orbiting an 8 Mq black hole, or a 0.5 
donor star orbiting a 1.3 Mq white dwarf — both of which would yield an inclination angle 
of ~ 60°. With regard to point (ii), it is possible that some evolutionary scenario other than 
the one utilized in our Monte Carlo study can explain persistent supersoft X-ray sources. 
We note that van Teeseling & King (1998) have recently proposed a model wherein the 
mass transfer and binary evolution in some supersoft X-ray sources are driven by wind loss 
from an X-ray irradiated, low-mass donor star. 

One difficulty with a small inclination angle arises in the effort to explain the 
photometric light curve ( |Will fc Barwig 1996| ). This curve shows a highly repeatable 
modulation at the binary period with a full depth of 0.5 mag. It seems difficult to 
explain such behavior with inclination angles as small as 20° — 30° (see, e.g., |Schandl 



[Meyer-Hofmeister fc Meyer 1997]) . Alternatively, the results could signify substantial 



vertical structure in the accretion disks of supersoft X-ray sources, perhaps as a consequence 
of the high rate of mass transfer. Some hydrodynamic disk simulations do show impressive 
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structures with significant vertical extent ([Blondin 1997 ). 

4.2. Kinematic Model for the Jet Lines 

We now examine more closely at least one possible interpretation of the jet lines. In 
this study, we interpret these lines as representing material flowing away from the accretion 
disk. In order to constrain the geometry of the jet, we adopt a very simple model for the 
outflow, which is based on spectral results (1) and (2) in §4.0. We define a jet cone into 
which the emitting atoms are released, all with speed Vj. The cone is defined by a tilt 
angle, 9j (with respect to the angular momentum axis of the binary), an opening angle 
of half- width, a, and an azimuthal angle, 0j. Note that 9j and 0j refer to the symmetry 
axis of the jet. Material is assumed to come off uniformly (per unit solid angle) within 
the cone during the times when the outflow lines are visible. There are jet cones located 
on either side of the accretion disk. In principle, the four jet parameters need not have 
identical values above and below the disk, and ideally we should flt for these parameters 
independently. However, since the blue-shifted line is unfortunately contaminated by P 
Cygni absorption, we focus our attention on the red-shifted portion of the outflow. 

At least flve effects contribute to the observed proflles of the jet lines. These include 
the Doppler shifts projected along the fine of sight due to (1) the velocity of the outflowing 
matter along the symmetry axis of the jet cone; (2) the orbital motion of the whole 
accretion disk around the center of mass in the binary system; (3) the different directions 
of atoms in the jet with respect to the axis of the cone; (4) the Keplerian motion of matter 
in the accretion disk around the white dwarf; and (5) the height above the disk of the jet 
emission region. The flrst two of these tend to dominate the position of the center of the 
line, while the next two tend to broaden the line and perhaps induce asymmetries. The 
flfth factor influences the detectability of binary motion in the velocities of the jet lines. 
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We first examine how the centroids of the fines vary witfi time around tfie orbit under tfie 
assumption that the peak of the fine is governed by (1) and (2) only. We later discuss 
calculations of the full line profile which justify this assumption, and we indicate what 
limits can be placed on effects (3), (4), and (5). 

With no loss of generality, we consider the case in which the binary orbit lies in the xy 
plane with the observer located in the xz plane at an inclination angle, i, with respect to 

— * — * — * 

the z axis. Let fc, V^, and Vj be, respectively, the unit vector pointing from the observer 
to the binary, the orbital velocity of the white dwarf, and the velocity of an atom moving 
along the jet axis with speed Vj in a frame orbiting with the white dwarf. In terms of their 
Cartesian coordinates, these vectors can be written as: 



where u is the Keplerian angular frequency of the binary, and 6j and (j)j are the instantaneous 
spherical polar angles describing the orientatiuon of the jet axis. With these definitions, 
and additional consideration of the systemic radial velocity (7) of the binary system, one 
can predict the observed Doppler velocity for the line centroids due to effects (1) and (2) 
above. If we define Vc and Vout as the Doppler velocities for the central and jet outfiow 
lines, respectively, we find: 



= [ — sin(i), 0, — cos(i)], 
Vo = Vo[ — cosiut), — sin(cjt), 0], 
Vj — Vj [sin(6'j) cos(0j), sin(6'j) sin(0j), cos(^j) ] , 



(4) 



(5) 



(6) 



Vc — ^ -\- Vo sin(i) cos(a;^) = 7 + cos(a;i) , 



(7) 



Vout = 7 + Kcos{ujt) — Vj sin(i) sm.{9j) cos(0j) — Vj cos(i) cos{9j). 



(8) 



We can now utilize our radial velocity measurements for RXJ0019 to extract the 
parameters associated with the outfiow. In order to proceed, we consider two possibilities: 
{case a) the jet axis is tilted by an angle 9j with respect to the normal to the accretion disk 
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(which is assumed to he in the orbital plane), and its projection into the orbital plane is 
fixed with respect to the line joining the two stars; and {case h) the jet axis is normal to the 
accretion disk, which in turn is tilted with respect to the orbital plane by an angle 6j. In 
both cases, the angle 9j is taken to be a constant; however, in case (a) the azimuth angle, 
0j, varies with the orbital period (i.e., = ut + constant), while in case (b) (f)j will vary 
with the disk precession time which, in general, should be much longer than the orbital 
period. 

In case (a) we can extract a rather tight constraint on 9j as follows. As shown above, 
the velocities of the central Ha line are consistent with the best fit for the Hell lines, 
i.e. Vc = VhbII- Therefore, we may effectively remove the first two terms in equation (|^) 
by subtracting VhcII from the observed Doppler velocities for the red jet line of Ha, i.e., 
Vout — VheII- We then compute the average of the subtracted results, which we associate 
with the last term in equation (P), i.e.: 

{Vout - Vueii) = Vj cos(i) cos{dj) = 815 km s"^ (9) 

Finally, we isolate the third term in equation (H) by first subtracting both Vneii and 
the result of equation (|^) from equation (H). We then fit a sine curve to the residuals of this 
operation. The result is: 

Vj sm{i) sm{9j) = 1.8±14 km [case (a)], (10) 

for the red outfiow line, where the fit was made to the data already folded at the orbital 
period. 

From the above analysis, we can set a combined limit on the inclination angle and the 
jet tilt angle by dividing equation (|10D by (||) to obtain, at 90% confidence: 

tan(i) tan(^j) < 0.027 [case (a)], (11) 
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for the red jet. If we now utilize the lower limit of ~ 16° for the orbital inclination angle, 
we can set a limit (again at 90% confidence) on the tilt angle: 

Oj < 5°.4 [case (a)]. (12) 

Returning now to case (b), where the azimuthal angle is assumed to vary only on 
a disk precession timescale, we analyze the radial velocities of the red outflow line in an 
analogous way to that just described. In this case, we assume only that the putative disk 
precession period is less than ~300 days. (For substantially longer precession periods, the 
only constraint that we can set on 9j, (pj, vj, and i is that the sum of the last two terms 
in equation (H) equals 815 km s~^ during the epoch of our observations.) For precession 
periods less than ~300 days (i.e., where (pj has moved through an angle much greater than 
27T over the interval of our observations), we can still isolate the middle term in equation 
(§) to set constraints on 9j and i. Again, we subtracted VhcII from the time series of red 
outflow velocities. We then carried out a Stellingwerf period search on the residuals of this 
operation. This search covered possible disk precession periods ranging from 0.5 to 300 d. 
From this analysis we derive a limit on the central term in equation (P) of: 

Vj sm{i) sm{9j) < 30 km s~"^ [case (b)], (13) 

at the 90% confidence level. In this case, the limit analogous to that given by equation (|TT]) 
is 

tan(i)tan(^j) < 0.037 [case (b)]. (14) 

Again, if we adopt a lower limit on the inclination angle of 16°, we find a corresponding 
limit on the tilt angle of: 

9j < 7.3° [case (b)]. (15) 

Finally, as noted previously, the blue outflow line at Ha is blended with P Cygni 
absorption, and the deconvolution of these features is beset with systematic problems 
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related to our limited spectral resolution and the uncertainties in the specific line profiles 
for each component. Nevertheless if we do investigate the radial velocity measurements of 
the blue outfiow line for evidence of a tilt angle with respect to the binary plane, as we have 
done for the red line, then the corresponding 90% confidence limits to equations (0) and 
(p!5|) are 10° and 15°, respectively, consistent with the results for the red line. 

From the above analyses, we conclude simply that the jet axes are essentially 
perpendicular to the orbital plane of the binary. 

4.3. Emission-line Profiles 

In the context of our simple kinematic model, we now consider two effects which would 
broaden the emission lines associated with the outfiow. Briefiy, these are (1) the finite 
opening (half) angle, a, of the jet, and (2) Keplerian motion of matter in the accretion disk 
that may feed the jet. We consider each of these in turn, and we assume in the calculations 
below that the outfiow emission region is optically thin. 

The opening angle of the jet affects both the width and the shape of the outflow lines. 
In order to constrain this opening angle, we adopt a simple model wherein the jet cone has 
a uniform flux of emitting atoms per unit solid angle over the entire cone. We then utilize 
equations (^) through to integrate the contributions over atoms moving in different 
directions within the outflow. However, equation (|^) must be generalized to refer to the 
velocity of any atom within the outflow, and we apply these considerations to a cone of 
atoms moving with a component along the red jet axis. Since we concluded in the preceding 
section that the jet cone axis is essentially normal to the orbital plane, i.e., 9j = 0, equation 
(I) becomes: 

14 = 7 + K sin{ujt) — Vj sin(i) sin(6') cos(0) + vj cos(i) cos(6'), (16) 
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where is the Doppler velocity of a single atom in the jet whose direction is given by 6 
with respect to the red jet axis. In order to compute the line profile due to the range of 
angles within the jet cone, one simply integrates the contributions from all atoms with 
6 < a (the half width of the opening angle), and < < 27r. After Doppler correcting the 
spectra to the rest frame of the accretion disk, one needs to sum only the contributions 
from the last two terms in equation ( ]TB| ) to find the line profile due to the finite opening 
angle of the jet. 

There is potentially a wealth of information in the profiles of the jet lines, which we 
earlier determined to have an intrinsic FWHM of ~ 400 km s~^. In principle, the central 
wavelength, width, asjTiimetry, and movement of the jet lines allow one to determine the 
outflow velocity as well as the origin and opening angle of the jet. However, given our 
modest instrumental resolution, we may only derive coarse constraints on these parameters. 

First, we calculated a grid of theoretical model profiles for various values of the opening 
half-angle, a, and the orbital inclination, i. For each grid point the jet line profile was 
computed and then convolved with the instrumental resolution (cr = FWHM /2.355 = 40 
km s^^ for the FLWO data). Next, the model profiles were renormalized in fiux and shifted 
in velocity to best match the observed profile of the red jet line, using a chi-square fit over 
the range of 6568 — 6620 A. The red jet line was isolated by fitting the central Ha line to a 
Gaussian and then subtracting the best fit from the spectrum. The line profile is distinctly 
asymmetric, falling off more steeply toward longer wavelength, as do the majority of the 
theoretical model profiles. Finally, we constructed contours of constant (minimum) chi 
square values in the a-i plane. The results are shown in Figure ^ Because the model is so 
simplistic, and after comparing the fits to the observations, we judge a reduced chi square 
value below 5.0 to be a reasonable representation of the observations; this contour is shown 
with a darkened line in Figure ^. Our simple jet model prefers a binary inclination angle 
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in the range of 35° — 60° with a jet opening half-angle in the range of 40° ^ a ^ 15°. 

If these rough limits on the opening half-angle are correct, then the jet outflows in 
RXJ0019 are not well coUimated. Alternatively, the jets could be well coUimated, while the 
line widths are increased by other effects, including a range of velocities within the line 
formation region. One particular additional cause of line broadening is motion of the disk 
material orbiting the compact star. The Keplerian velocity at radial distance r is 



where rio is the radial distance from the compact object (assumed to be a white dwarf) in 
units of 10^° cm. If the accretion disk lies in the orbital plane and the system is viewed from 
an inclination angle, i, then the maximum Doppler velocity due to matter orbiting in the 



from the summary of spectral facts cited at the beginning of §4, we see that the FWHM of 
the Ha hues (central and outflow hues) is only 400 km s~^. Therefore, we conclude that the 
regions of the accretion disk contributing to the Ha emission are located relatively far from 
the the white dwarf at radial distances 



For a white dwarf of 1 Mq and a radius of ~5 x 10^ cm, the Ha emission region must lie 
at least ~40 white dwarf radii from the compact star. This suggests that the 815 km s~^ 
jet is not produced near the compact star, but rather emanates from a region further out 
in the disk. The escape velocity from a 1 Mq white dwarf is ~7000 km s^^, so that the 
observed projected outflow speed of 815 km s~^ again suggests that the ejection is taking 
place relatively far from the surface of the white dwarf. 

A final caveat in regard to possible broadening due to the Keplerian motion of the disk 
material is in order. If the outflow originates very near the compact star, e.g., if the outflow 
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accretion disk in an annulus at radius r is ±1150 {Myj^./ Mq^^'^Tiq^'^ sin(i) km s 



. However, 




(18) 
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were jet-like, and the mechanism driving the outflow could remove the Keplerian motion of 
the inner disk from the ejected matter, then perhaps the observed Doppler velocities would 
not have any broadening due to the rotation of the disk material. In this case, our limit for 
the inner radius of the emission line region given by equation (|18D might not be valid. 



Lastly, we comment on the fact that both the red and blue outflow emission lines, 
when present, always have approximately the same intensity. This implies that the emission 
region must be far enough from the disk that the outflow on either side is equally visible. 
On the other hand, the emission region cannot be so far from the disk that the orbital 
Doppler effect that we clearly detect, would be smeared out by emission originating over 
different orbital phases. 

The flrst requirement implies the obvious geometric constraint on if, the height above 
the disk at which the emission takes place. We may express the constraint on H relative to 
the semimajor axis of the orbit of the compact star, a^.: 

»>^^^R, (1 + ,) aj6£(l+^ 

ax axtaji[i) aq tan(z) tan(«) q 

where Rd is the radius of the accretion disk, q is the mass ratio Mdon/Mwd, and / is 
the fraction of the Roche lobe of the compact star that the accretion disk occupies. The 
approximation for Rn/a in equation (ITW) is accurate to within ~ 10% over the range q < 2. 



For typical system parameters and i ~ 25° and / = 0.7, we expect H/ax ^ 1. 

The Doppler modulation of the jet lines at the binary period constrains the upper limit 

for the emitting region above the accretion disk. This limit can be expressed as: 

/ 815 km s^^ 
H < v,/uj = . (20) 

We can also cast this limit in dimensionless form by dividing by ax, the semimajor axis of 
the orbit of the compact star: 

H Vj Vi , ^ , , 

- < = TTI = 11-5 tan z . 21 

ax ooax K/smyi) 



-27- 



For i ~ 25°, we find that H/a^ < 5.3. 

In summary, if the emitting region is no further above or below the accretion disk 
than ~6 ax, nor any closer than ~1 a^, then both the red and the blue jet lines should be 
approximately equally visible while retaining the net Doppler motion at the binary period, 
as we have observed. 

5. P Cygni Profiles 

As shown in Figure |^ the strength of the P Cygni absorption varies with orbital phase. 
The strength of the absorption feature is a minimum at spectroscopic phases 0.8 to 1.0, 
when the central line is moving from its rest wavelength position to its maximum redshift. 
A similar result was reported by [Reinsch et al. 1996b| . If the shape of the P Cygni trough 
and its central wavelength are independent of orbital phase, then the Doppler motion in the 
central line could reduce the net absorption strength near phase 0.5. Since that is nearly 
the opposite of what we observe, we conclude that the phase dependence of the absorption 
strength is intrinsic to the geometry of the absorbing gas. 

The nature and origin of the gas responsible for the P Cygni absorption remain a 
mystery. Our first hypothesis was that the gas responsible for the jet outfiow might also 
be the gas that causes P Cygni absorption. The phase dependence of P Cygni absorption 
implies an opening half-angle ~ 100° for the absorbing gas which is much larger than the 
~ 40° upper limit for a in the jet. Furthermore, the absorption cone is not only wide, but it 
is azimuthally asymmetric, while the jet outfiow has been shown to be closely aligned with 
the binary spin axis. Finally, the jet and P Cygni components appear to have differences 
in their persistence time scales. Although the P Cygni absorption strength is variable, it is 
apparent in all of the observing runs. By contrast, the jet lines are present only in about 
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one-third of our observing epochs. Therefore, the very different behavior between the jet 
hnes and the P Cygni absorption features leads us to beheve that they arise from different 
physical regions in the binary system. 

6. Summary and Conclusions 

1. We have made spectroscopic observations of RX J0019. 8+2156 over a 2.5 year 
baseline, collecting over 200 ks of data on this source. 

2. Radial velocity studies of the Hell A4686A emission line have yielded a Doppler curve 
with a velocity half-amplitude of 71.2 ± 3.6 km s~^(see Table The binary mass 
function is 0.0247 ± 0.0040 M©. 

3. The period as determined from our radial velocity study is in agreement with the 
photometrically determined value obtained from the literature. However, the phases 
determined by the two methods disagree by 28° {6a significance), assuming that the 
binary orbit is circular and the accretion disk is circularly symmetric. 

4. The summed rest-frame spectrum shows extraordinary Hell series emission. Most 
notably, the n— i>5 transition series is visible for values of n up to 22. 

5. A model-dependent population synthesis study indicates that if the binary parameters 
are to be consistent with the measured mass function, then 16° ^ i ^ 40°. 
Furthermore, if the system also has a mass ratio which is capable of driving high 
mass transfer rates (in the context of the van den Heuvel et al. model), then the 
orbital inclination angle is 16° <^ i ^ 30°. However, such a low inclination is difficult 
to reconcile with published photometry, which shows a pronounced primary and 
secondary dip. 
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6. We discovered the existence of low-velocity bipolar jets that are revealed as red-shifted 
and blue-shifted emission lines that flank Ha, H/?, and Hell A4686. These lines are 
transient and were detected only during the 1994 Dec, 1996 May, and 1996 Dec. After 
correcting for the gamma velocity of the system, these lines are nearly symmetrical 
around the central line with a velocity of ± 815 km s~^. 

7. The outflow lines also exhibit orbital Doppler motion with the same amplitude and 
phase as that determined by the HeH lines. This imphes that the jets are oriented 
nearly perpendicular to the orbital plane. 

8. The profiles of the outflow lines potentially possess a great deal of information. 
While the limited spectral resolution limits our analysis, we are nevertheless able to 
conclude that the range of inclination angles is hkely to be 35° ^ i ^ 60°, and the 
corresponding range of opening half-angles 40° ^ a ^ 15°, if the intrinsic width in the 
line profiles is dominated by the opening geometry of the jet. 

9. In order for both red and blue outflow lines to be observable at all orbital phases and 
also exhibit the Doppler curve of the white dwarf, the emission must originate from a 
distance H from the disk in the range 1 ^ H/ux ^ 6, where is the semimajor axis 
of the orbit of the compact object. 

10. We discovered that the strength of the P Cygni profiles is variable with orbital phase. 
The ~ 100° opening half-angle of the absorbing gas and its asymmetry with respect to 
the line joining the stars is distinctly different from the geometry of the gas emanating 
from the jets. For this reason, and because of the transient nature of the outfiow lines, 
it is unlikely that the absorption and jet components have the same origin. 
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Figure Legends 

Fig. 1. — Stellingwerf period search using radial velocity measurements of RXJ0019. Minima 
in the variance statistic correspond to periodicities in the data. 

Fig. 2. — Enlargement and fit to the minimum in the variance statistic near a trial period 
of 15.85 hr. 

Fig. 3. — Hehocentric velocities of the Hell line folded about the 15.85151 hr period and a 
sine wave fit. The data span an interval of ~2.5 years. 

Fig. 4. — Top: Deep spectrum of RXJ0019 in the emission-line rest frame (see text), prior to 
continuum removal. Bottom: Enlargement of the 5600-6450A portion of the deep spectrum 
after continuum subtraction. Most of the emission lines are the n— >5 transitions of He II (see 
Table 4). 

Fig. 5. — A fifteen minute observation (per red and blue exposure) of RXJ0019 made with the 
FAST spectrograph and the FLWO 1.5°^ telescope on 12 December 1994. Signatures of high 
velocity lines and P Cygni type absorption are readily seen. The insets show enlargements 
of the Hell A4686, H/3, and Ha profiles (left to right). 

Fig. 6. — Velocity of the red outflow (top), Ha (middle), and blue outflow (bottom) emission 
as a function of the spectroscopic phase. The Ha line is used because it has the strongest 
outflow components. The solid lines show the best fit sine curves whose parameters are 
summarized in Table 5. 

Fig. 7. — Equivalent width of the H/3 emission (bottom), P Cygni absorption (middle), 
and the ratio of the two ( top ^ as a function of the spectroscopic phase. The dots (triangles) 
represent observations with (without) observed high velocity lines. H/? lines are used because 
the absorption trough is better defined. 
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Fig. 8. — Histogram of the number of simulated systems as a function of final orbital 
inclination (left panels) and final donor-star mass (right panels). For each case we separate 
the systems that reach a 15.85 hr binary period early in their evolution (P < 0, top) from 
those systems that have evolved past their minimum period (P > 0, bottom). The unshaded 
histograms represent systems with mass transfer rates below 3xl0~^ Mq yr~^, while the two 
levels of shaded histograms represent mass transfer rates above 3x10^*^ and 10^^ Mq yr~^ 
respectively. 

Fig. 9. — Contours of the best-fit value of reduced chi square as a function of opening 
half-angle of the jet and the inclination angle of the binary system. 
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Table 1. Observations of RXJ0019.8+2156 at MDM 



Observing JD-2400000 JD-2400000 Exposure Bandpass Number of Number of Number of 
Run Start End (sec) (A) Nights Pointings Spectra 



Dec 94 


49696.693 


49700.645 


300 


3800-6200 


4 


4 


33 


Jun 95 


49875.973 


49881.932 


100 


4100-7200 


6 


8 


88 


Oct 95 


49991.667 


49996.958 


200 


4600 7575 


6 


28 


140 


Nov 95 


50025.694 


50026.711 


200 


3885 6300 


2 


4 


23 


Dec 95a 


50050.580 


50051.802 


300 


4500-7425 


2 


6 


28 


Dec 95b 


50052.561 


50055.740 


300 


4025-7110 


4 


13 


60 


Dec 95c 


50056.646 


50065.738 


300 


4280-7125 


10 


18 


115 


Jan 96 


50086.572 


50094.566 


300 


4280-7140 


9 


14 


87 


May 96a 


50236.985 


50241.967 


200 


4521-7485 


6 


6 


17 


May 96b 


50249.974 


50270.954 


200 


4536-7485 


22 


16 


91 


Dec 96 


50448.601 


50454.607 


300 


4285-7150 


6 


6 


48 


Feb 97 


50465.630 


50470.566 


200 


4460-7400 


6 


6 


27 
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Table 2. Observations of RXJ0019.8+2156 at FLWO 



JD - 2400000 Exposure Bandpass 
(days) (sec) (A) 



49655.738160 


1200 


4569- 


-5520 


49656.755360 


900 


4574- 


-5344 


49656.792080 


900 


5733-6687 


49662.797710 


900 


4570- 


-5523 


49662.819210 


900 


5729-6683 


49686.570080 


900 


4572- 


-5523 


49686.588220 


900 


5738- 


-6684 


49689.720000 


900 


4571- 


-5522 


49689.748870 


900 


5733- 


-6687 


49694.689250 


900 


4572- 


-5523 


49694.709460 


900 


5734- 


-6688 


49695.694650 


900 


4570- 


-5521 


49695.715140 


900 


5733- 


-6686 


49720.573440 


900 


4574- 


-5529 


49720.616370 


900 


4574- 


-5528 


49720.624670 


300 


5741- 


-6687 


49720.653580 


900 


5741- 


-6687 


49724.613470 


900 


4573- 


-5524 


49724.633840 


900 


5736- 


-6689 


49747.577730 


900 


4571- 


-5522 


49747.597220 


900 


5694- 


-6682 


49751.597410 


900 


4572- 


-5523 


49751.617640 


900 


5752- 


-6727 


49754.580720 


900 


4573- 


-5528 


49754.594210 


900 


5744- 


-6725 
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Table 3. Orbital Parameters for RXJ0019.8+2156 



Parameter 


Results 


Reference 


P (Hell velocities, days) 


0.6604795 ± 0.0000096 


This work 


Ts (max. velocity, HJD) 


2450002.2713 ± 0.0046 


This work 


7 (systemic velocity, km s~^) 


-57.5 ± 2.3 


This work 


K (velocity semi-amplitude, km s~^) 


71.2 ± 3.6 


This work 


Mass function (M©) 


0.0247 ± 0.0040 


This work 


P (photometric, days) 


0.6604721 ± 0.0000072 


Will & Earwig 1996 


To (photometric min., HDJ) 


2450014.273 ±0.007 


Will & Earwig 1996 


P (photographic, days) 


0.6604565 ± 0.0000015 


Greiner & Wenzel 1995 


To (photometric min., HDJ) 


2435799.247 ± 0.033 


Greiner & Wenzel 1995 
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Table 4. Lines identified in the rest frame spectrum of RXJ0019. 



Obs A Species Rest A EW 



3924, 


.8 


He II (15-4) 


3923.5 


0.29 


3976, 


.9 


He 


3970.1 


0.22 


4026, 


.3 


Hell (13-4) 


4025.6 


0.62 


4109, 


.3 




4101.7 




4199, 


.7 


Hell (11-4) 


4199.8 


0.68 


4343, 


.5 


H7 


4340.5 


0.43 


4542, 


.1 


Hell (9-4) 


4541.6 


0.94 


4686, 


.1 


Hen (4-3) 


4685.7 


9.72 


4862, 


.3 


H/3 


4861.3 


1.61 


4930, 


.4 






0.33 


4947, 


.7 






0.15 


5290, 


,2 


OVI 


5291 


0.90 


5412, 


.0 


HeH (7-4) 


5412.0 


1.74 


5665, 


.2 






0.14 


5803, 


.0 


CIV 


5801.5 


0.31 


5814, 


.2 


CIV 


5812.1 


0.12 


6004, 


.5 


HeH (22-5) 


6007 


0.11 


6037, 


.2 


HeH (21-5) 


6036.7 


0.11 


6074, 


.3 


HeH (20-5) 


6074.1 


0.17 


6117, 


.8 


TT TT /1 n cr \ 

Hell (19-5) 


6118.2 


0.23 


6170, 


.6 


HeH (18-5) 


6170.6 


0.18 


6200, 


.7 






0.24 


6234, 


.3 


HeH (17-5) 


6233.8 


0.32 


6311, 


.7 


HeH (16-5) 


6310.8 


0.30 


6374, 


.5 


FeX 


6374.5 


0.52 


6406, 


.9 


HeH (15-5) 


6406.3 


0.43 


6502, 


.7 






0.42 


6528, 


,1 


Hell 


6527.9 


0.60 


6543, 


.9 






1.1 


6563, 


.2 


Ha 


6562.8 


13.84 




o 






1 



Notes 



Almost all emission absorbed in P Cygni 



Almost all emission absorbed in P Cygni 



Strong P Cygni Profile 



Strong P Cygni Profile 



Part of blend? 



At edge of absorption band 



P Cygni absorption present 
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Table 5. Fitted parameters to the emission-line velocity curves. 



Feature 7^ 7 — ^Heii^ Phase^ K"^ 





(km s 




(km s 










(km s" 




Red 


758.0 ± 


10.0 


815.5 ± 


10.0 


0.787 


± 


0.032 


73.1 ± 


14.5 


Central 


-40.7 ± 


7.2 


16.8 ± 


7.2 


0.809 


± 


0.033 


48.9 ± 


10.44 


Blue 


-909.2 ± 


14.9 


-851.9 ± 


14.9 


0.775 


± 


0.032 


105.5 ± 


20.7 


Hell A4686 


-57.5 ± 


2.3 


0.00 




0.750 


± 


0.007 


71.2 ± 


3.6 



'^The quoted error bars represent the la confidence limit after scaling the error 
bars on individual measurements such that — 1. 
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09 05- OOL- OQL- 

(s/lu>]) A^ioo|a/\ oij^u930i|aH 
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